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Summary 


Unmanned aerial vehicles (UAVs) are critical components of the future naval forces. UAV 
control and monitoring with autonomous operation will become an absolute necessity and 
adaptive cooperation of vehicles is the only practical alternative. 

The objective of this project is to develop and evaluate new methodologies for cooperative 
(formation) control of multiple unmanned air vehicles. The goal is to have multiple UAVs 
working together as a group. Instead of separately assigning distinct tasks to each vehicle, the 
operator would assign tasks to the UAV group, which then determines the best way to 
accomplish each task, freeing the operator to maintain surveillance over the entire operation. In 
this project we investigated Path Tracking and obstacle avoidance of UAVs using fuzzy logic 
method. Algorithms for close formation control of multi-UAVs are developed and simulated. We 
also investigated fault-tolerant control of single UAVs by neuro-adaptive method. Detailed 
description of this method is provided in this document. The project has supported 5 graduate 
students with 9 technical papers published. 

Technical Details 

A. Path Tracking and Obstacle Avoidance of UAVs 

Unmanned Aerial Vehicles (UAVs) have attracted increasing attention in the military and 
civilian applications, such as mapping, patrolling, search and rescue, and reconnaissance. In 
these types of applications, the Path tracking and Obstacle avoidance are fundamental. In this 
work, we developed a fuzzy logic based approach to path tracking and obstacle avoiding. We 
consider the case that obstacles are either still or moving and appear along the pre-determined 
flight path unexpectedly. By using suitable sensors we identify the relative distance between the 
vehicle and the obstacle and make timely adjustment to the pre-planned flight path. Fuzzy logic 
control algorithms are developed to achieve close path tracking while avoiding obstacles. Some 
simulation results are presented in Figure 1-Figure 6. 



Figure 1 - Path tracking with one obstacle - shape 1 
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Figure 3 - Path tracking with one obstacle- shape 3 



Figure 4 - Path tracking with two obstacles along the same Line 
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Figure 5 - Path tracking with two obstacles across the flight path 


B. Fuzzy Logic based Control of multi-UA Vs 

We also explored a method to achieve close formation tracking control of multi-UAVs by 
applying fuzzy logic theory. The formation trajectory is achieved by the control of Wingman’s 
heading velocity and heading angle while the lead UAV is maneuvering the given mission path. 
Since the values for specifying linguistic terms of the controller output for trailing UAV actually 
determine the shape of the input-output shape, they are selected as the main parameters for 
optimization. We have successfully used such a fuzzy controller to approximate the nonlinear 
system, which can be applied for the supervised learning. 

The method guarantees the desired tracking precision and also ensures the control action 
smooth under system uncertainties and external disturbances. Simulation of leader-wingman 
formation geometry was conducted. Both theoretical studies and simulation results demonstrate 
that this control algorithm is simple, less sensitive to the flight aerodynamic model. The 
formation geometry is determined by the relative position between the Leader and Wingman as 
shown in Figure 6. The formation control objective is to steer the Wingman (follower) to 
maintain certain separation distance in longitudinal, lateral and vertical directions. 



Figure 6 - Formation flight in inertial frame 
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Table 1: Parameters Used for Simulation 


Wing Area 

300 

fa 

Wing Span 

30 

ft 

Mass 

776.4 

lb 

Altitude 

45,000 

ft 

Dynamic Pressure 

155.8 

lb/ft2 

Mach 

0.85 

mach 

[x, y, z] 

[lOOft 60ft Oft 
[ 90ft -50ft 0ft 


[x d ,y d , z d ] 

'60ft 30*/4ft Off 
60ft -30^/4 ft 0ft 


[V w (0\ <M0X MO)] 

[750ft/s, n/6 , 30,000fl] 

[V L , n. A1 

[825ft/s, */8, 30,000ft] 



Figure 7a - Tracking trajectory of a sinusoid function 




Figure 7b - Separation tracking errors of wingman A 
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Figure 7c - Separation control signals of wingman A 
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Figure 7d - Separation tracking errors of wingman B 
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Figure 7e - Separation control signals of wingman B 
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Neuro-robust and adaptive control algorithms and memory-based control algorithms for close 
formation control have been developed. A series of simulations have been carried out. Some of 
the results based on the parameters and conditions as in Table 1 are presented in Figure 7a- 
Figure 7e. As clearly indicated from the results, the proposed method maintains a desired 
tracking precision and also ensures the control action smooth under the system uncertainties and 
external disturbances. 

C. Neuro-adaptive based Fault-tolerant Control 

One of the typical yet devastating faults in flight systems is the loss of the driving (propulsion) 
power. In this work, we explore an innovative approach to accommodate such faults in 
unmanned aerial vehicles operating under varying fight conditions. 

We present a neuro-adaptive control strategy, composed of two Neural Network (NN) units, 
the former to adapt to unanticipated uncertainties and the latter to enhance performance by 
compensating (funnel ing) the estimation error caused by the first NN to incorporate the 
important and most common problem faced by the airborne vehicles, unforeseen and uncertain 
actuating failures. New algorithms are derived to cope with the sub-system failures due to the jet 
engine partially losing its propulsion power ensuring the system stability. Simulation study using 
a generic model demonstrates a dramatically improved performance in the face of fading power 
faults and system uncertainties. Some of the details are presented as follows. 


C.l. Introduction 

Every component in a flight vehicle provides a certain vital function and the overall system 
works satisfactorily only if all components provide the service they are designed for, accurately. 
Therefore, a fault in a single component usually degrades the performance of the overall flight. 
In order to avoid performance deteriorations or damage to the vehicle, faults have to be identified 
and accommodated in a timely manner. In UAV, failures of control elements (sensors, actuators, 
etc.) are often encountered. A weak element in this framework is control loops. In fact automated 
systems are vulnerable to faults such as defects in sensors or in actuators, which can cause 
undesired reactions and consequences. A conventional feedback control design for a complex 
system such as UAV may result in unsatisfactory performance, or instability or malfunctions in 
system components, causing serious operation and safety problems. Therefore, fault tolerance 
has been one of the major issues in control of flight vehicles. 1 . 

A closed-loop control system able to tolerate component malfunctions while maintaining 
desirable system performance and stability exhibits fault-tolerant property. These systems, also 
termed as self-repairing systems, have been a subject of widely scattered studies. In this 
framework, the design of Fault Tolerant Control architecture is of crucial importance and 
solutions aiming at adapting the control strategy to the presence of the fault are needed in order 
to achieve prescribed performances. This is usually achieved by providing the control loop with 
a decision making layer that analyzes the behavior of the plant and adapts the control strategy to 
hold the controlled system in a region of acceptable performances. This is also phrased as Fault 
Tolerant Control (FTC). 

Fading power failure causes major problems in many critical control systems such as flight 
control systems. It is often not known when a jet propulsion fails and how much the failure is but 
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the remaining power could be still enough to accomplish a desired control task such as 
emergency landing of a flight vehicle. The question is whether a control system is intelligent 
enough to use remaining thrust in the presence of unknown failures. Figure 8 illustrates the 
general fault scenario. An adaptive approach, which is capable of controlling systems with 
uncertainties, is thus of interest in developing control schemes which are efficient for handling 
unknown sub-system failures. 

Most existing fault-tolerant control systems are built with redundant controllers 24 . The 
backup control loop is employed once failure of the main controller is detected online. However, 
often times, the additional cost, space and complexity of incorporating redundant hardware 
makes this approach unattractive. The use of redundant controllers may significantly increase the 
capital and maintenance costs of such control systems to unacceptable levels. Therefore, fault- 
tolerant control systems with no or less control redundancy are desirable. 


Redundancy free fault tolerant control has been studied recently. One approach is to detect 
and diagnose system faults online and update control laws accordingly to stabilize faulty 
system 1 . However this method is usually difficult to implement, as it required multiple control 
algorithms to be pre-programmed for each possible failure scenario. Various approaches to 
achieve control reconfiguration have been researched 5, 6 ' 8 . Fault-tolerant control systems can 
benefit from blending the perspectives of intelligent control. Passin o 9 describes some of the 
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Figure 8 - The general fault scenario 


limitations of current approaches, and states that a combination of intelligent and conventional 
control methods may be the best way to implement autonomous control. However, from the 
application point of view, only a few studies have been devoted to flight vehicles. Rauch 10 talks 
about fault accommodation in flight and space shuttle control, where Probabilistic Neural 
Networks are used. Polycarpou et al. 11 presented a neural network (NN) based fault tolerant 
approach and application on Beaver aircraft but did not mention fault occurrence in the primary 
flight control system. Chen et al. 12 presented an LMI-based synthesis scheme for robust fault- 
tolerant control systems design, where an unmanned aircraft suffering from wing impairment is 
used as a design example. Boskovik et al. 13 presented a new Decentralized Failure Detection and 
Identification (FDI) and Adaptive Reconfigurable Control (ARC) scheme, well suited for 
achieving the desired flight performance in the presence of multiple and applications to F-18 
aircraft carrier landing maneuver. NN model based FTC has been studied in recent years 14 . 
However, the key point of applying NNs to the FTC of a practical system is that this model must 
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be able to approximate the system faults with sufficient accuracy and satisfy the real time 
requirement. 

Our objective is to design a fault tolerant flight control system guaranteeing stability and 
satisfactory performance not only when all components are operational, but also in the case when 
actuating power fades and failures occur in flight dynamics. We propose a neuro-adaptive 
control scheme that is capable of tolerating component malfunctions and fading jet-engine 
power, while still maintaining desirable performance and stability properties. It possesses loop 
integrity, reliability, maintainability and survivability. The NN used in this work is an online- 
trained network in which the weights are adjusted online accordingly to accommodate faults 
occurring during system operation. 

C.2. Flight Vehicles Dynamics with Faults 

FTC systems can be characterized as robust, reconfigurable or a combination of both. The flight 
vehicles, regardless whether they are military or commercial, require precise maneuverability 
and control; therefore, the design of a robust controller that assures a safe and precise operation 
of the flight vehicle is of extreme importance. A great difficulty related with the design of any 
flight control system is that a complete and accurate dynamic model - including aerodynamic 
coefficients, is required. However, it is difficult to identify accurately the aerodynamic 
coefficients because they are nonlinear functions of several physical variables. Neural networks 
have been proposed recently as an adaptive controller for nonlinear systems. The development 
and analysis of NN based control algorithms for the control planes of flight vehicles with faulty 
dynamics, varying mass, nonlinear velocity and parameter uncertainties is a very challenging 
task. 

The trajectory tracking of a flight vehicle is achieved through the help of control planes, 
which in turn are adjusting the angles of attack, sideslip and body-axis roll of the flight vehicle. 
These angles can be observed in Figure 9. We use the dynamics similar to X-37 represent the 
nonlinear mathematical model of flight vehicles. The controllers are designed with the 
assumption that some of the aerodynamic coefficients are fully understood. 



Figure 9 - Coordinate System 
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The body-fixed axes, nonlinear equation of motion for a UAV over a flat Earth are given by 
Lee and that describes a similar model as follows 15,16 : 
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Definition of state and control variables, forces and moments in the preceding equations are 
described in the Nomenclature. It is assumed that the aerodynamic forces and moments are 
expressed as functions of angle of attack, sideslip angle, angular rates, and control surface 
deflection. L, M, and N are expressed as follows 16 : 


L = 

M = 

N = 
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The moments, coefficients and parameters are described in the Appendix. Substituting the 
aerodynamic coefficients into the flight dynamic equations yields the following: 
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Naming the states x = [a p ^ , y = [p q r] r and z = [9 , the nonlinear dynamic 

system under consideration can be described as follows: 

( 10 ) 


x = A(a,PJ,0)+ B(a,p,<f>,9)y{p,q,r) 
y = C(a,P,p,q,r) + D(a)y(p,q,r) + E(a,P)u 
z = H(t,d)y(p,q,r ) 


( 11 ) 

( 12 ) 
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where A, B, C, D, E and H represents the grouped terms of Eq. (10) - Eq. (12) in order and 
«=[*. s, s r ] is the control input 


Remark: 

The mass m of the reusable vehicle is considered time-varying due to the m f fuel 
consumption - especially during the ascend mode. Therefore, the total mass of the vehicle is: 


m = m 0 + m f + m p 

where 

m 0 real dry weight 

m f fuel weight (variable) 

m p payload weight 

Note that this change in mass affects the dynamics, however, the change in UAV’s mass dm/di is 
considered unknown, making m unknown in the model, thus should not be used directly in 
control design. 

Observations: 

• The simulation studies show that the norm of the gain matrix from Eq. (10) is negligible, 
and therefore, this leads to a simplification of the dynamic model under consideration 
without any lost of generality. 

• The dynamic system equation described in Eq. (10) - Eq. (12) is controllable only if the 
matrix BE is invertible. Matrix E is known from the literature as being invertible. As for 
matrix B, the control surfaces of the reusable launch vehicle are designed to control each 
axes angular rate of aircraft independently; therefore, the matrix B is invertible for all 
cases. Numerical studies for the aerodynamic model also show that B is always invertible. 
The dynamic system presented can be recast in the second-order form by deriving Eq. (8) 
and substituting Eq. (9); it can follow: 


x — — A (.) + BC + 
dt W 


— B(.) + BD 
dt W 


y + BEu 


(13) 


The above equation is the governing equation of motion for the aircraft and will be used for 
the control design and stability analysis processes of the following sections. All the effects 
caused by possible failures in the system can be lumped into an uncertain term. The fault 
dynamics can be described as follows: 

f f = ±x(t-T,) ru with 

fo if t<T, (14) 

XU-T f ) = \ ’ 

1 - e ' ' otherwise 
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Note that the faulty function described above covers both the jump and continuously time 
varying faults in both directions. The additional term added to the system dynamics now 
represent the effect of system faults. The model of the actuator power loss failure can be 
represented as 17 

u f = E f u (15) 


E f e R'”’ is the actuator fault matrix with the following diagonal structure 

e /t 0 .... 0 

0 .... 0 

. 0 

0 0 0 


where 

= 0 represents power loss 
0 < s fi < 1 represents fading power 
e fi = 1 represents healthy actuator 
i = 1,2 


(16) 


Properties 

The faulty dynamics defined in Eq. (15) have the following properties. 

• The matrix £/is a diagonal matrix that is positive definite. This fact can be verified by 
observing the limits of e fi . 

• Such matrix is nonsingular; therefore all the eigen values are non-zero. 

• It is real symmetric; therefore all the eigen values are real numbers. 

• The eigen values of such matrix can be given as : 

E f =diag(s f i. s^) 

=> £ fV ...e fn are eigen values of E f 

Thus, the dynamics of the RLV Eq. (13) become 


Define error dynamics as 


x = — A(.)+BC + 
dt 


— B(.)+BD 
dt 


y+RE/U+f, 


e-x-x* and s = e + k 0 e 

where ko> 0 is a design constant. The control gain is represented as follows 

R=BE = R 1) +AR 


(17) 


(18) 
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Error dynamics of Eq. (17) can be simplified with the error function defined as in Eq. (18) and 
represented as: 

s = f 0 +W+RoE f u (19) 


where 


/ 0 = BC+BDy-x + k 0 e- Known, Time variant 
d d 

Af = — A(.)+ — B(.)y+ARE f u ± /, - Unknown, Timevairant 

dt dt (20) 

Rq - Known Part of control gain 
A R - Unknown Part of control gain 

Note for the sake of simplicity, the arguments of the functions were omitted, even though this 
does not affect the clarity. Also a part of control gain Ro is available if not the total control gain 
REf. The framework for a control to deal with this type of failures is illustrated in Figure 10. 


Actuator System 

Faults Faults 



Figure 10 - Controller framework for fault tolerance 

Assumptions: 

For the system to admit a feasible solution under faulty condition, we need to assume that 
there exists a positive constant c 0 > 0 such that ||A£|| + c 0 < ||/?|| 


Properties: 

1) Matrix E is known from the literature as being invertible. As for matrix B, the control 
surfaces of the UAV are designed to control angular rate along each axis of the aircraft 
independently; therefore, the matrix B is invertible for all cases. Numerical studies for 
the aerodynamic model also show that B is always invertible. 

2) £/' = E/ t > 0 

3) If Ro is full row rank, then there exist constants c mm and such that 

C mi. IHf ^ s T RoE f R/s < Cnax HI 2 
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where c m „ > c mm > 0 denotes the maximum and the minimum eigen values of RoEjRo 
respectively. This proposition is valid since the matrix is symmetric and positive-definite. 
This can be justified by noting that £/ is a diagonal matrix and Ro = BE is invertible as 
explained above. Also observing the properties of £/ defined in the previous section, the 
matrix has positive, real eigen values. 

Now we proceed to design an adaptive control algorithm for the above described highly 
coupled, nonlinear, uncertain, faulty complex system. 

C.3. Neuro-adaptive Fault-tolerant Control Design 

The major challenge in control design for the above described system is the nonlinear and 
uncertain nature of A R and Af The basis for most stable adaptive controllers has been to 
reorganize Af as a linear combination. However, this method is not suitable for complex systems 
such as aircrafts since it involves complicated design procedures and heavy computations. 
Another approach is to estimate Af via multilayer NNs and train these networks via gradient 
algorithms. Recognizing the potential instability associated with these algorithms, several 
researchers have proposed NN control schemes using on-line and stable training mechanism 
based on Lyapunov stability theory. We further extend these works in an innovative path. 

We “reconstruct” Af via the first NN unit as 

fm = /o + Vi 

Where 

fo the available part of f 

Wj e R mxm the optimal weight of the NN 

V, e R m the base function of NN whose selection typically include linear function, sigmoid 
function, hyperbolic tangent function, radial basis function, etc. 

Let the discrepancy between / and fm be denoted by: 

£ -f~ fm (21) 

where e is the reconstruction error. With the support of the approximate theory of multilayer 
neural networks, such an error can be made arbitrarily small provided that certain conditions 
(e.g., sufficiently large numbers of neurons, sufficient smoothness of/ etc.) are satisfied. This 
has been the basis for many NN based control laws. In reality, however, it should be noted that 
one can only use a network with a finite number of neurons and that the nonlinear function being 
reconstructed is not smooth enough. Thus, during system operations, the reconstruction error 
sometimes is bounded by a constant, but at other times it may not necessarily be confined to that 
number. For this reason, simply assuming that ||e|| is uniformly bounded by a constant does not 
warrant the success of the control scheme, as seen in Figure 11 . 


14 




Figure 11 - The possible profile of NN reconstruction error 


For executing the idea, in our development, we remodel reconstruction error as follows 

e = ce, + (1 - c)e 2 

\ 1 if the NN works satisfactorily 

with c — { 

[0 if the NN fails to work well 

To make the problem solvable, we make the assumption that 

ll^ll^cco, andKIMI/-/.-<ViMO. 

The point to be made is that do and d\{.) are still unavailable and should not be included in the 
control scheme directly. Another major challenge due to the fact that we do not know when the 
NN behaves well. Sanner and Slotine 19 proposed a method of partitioning the state space into Aj 
and a' and assumed that the NN works satisfactorily within Ad. Here we take a different method 
to address this issue. Namely, we incorporate two NN units in the control scheme, with the first 
NN unit compensating the lumped nonlinearities/uncertainties and the second NN unit 
attenuating the effects due to the NN reconstruction error and the other resulting uncertainties. 
The overall control scheme is given by 

u = R v k 0 s— (22) 

The estimate neural network of the uncertain partzf/is: 

/*-/. + *>, (23) 

The terms used in the controller and estimation of the first neural network are: 

u ca = a compensating signal to be specified 

ko> 0 = a constant chosen by the designer 

w, = is the estimate of the ideal weight values 

Before defining the signal u ca the following signals and scenarios are to be considered. It is 
interesting to observe that the control algorithm remains unaffected even though the control gain 
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is considerably affected by the actuator failures. This is one of the features of the neuro-adaptive 
controller. It adapts to the changes in the uncertainty without necessitating the changes in the 
control algorithm. With the control of the above scheme, the original system Eq. (19) becomes 

= -K^ E fK S + ^V, - K E fK U ca +7 (24) 

where W, = W ] -W ] is weight estimation error and 

7 = £ + (/ - R V E / Rj ) (/ 0 + W? V,) (25) 

is the resultant uncertainty to be attenuated by the second NN unit. Expressing 7 in a similar 
fashion as done earlier, we have, || 7 || < y a 

N,+ 3 

r. = L (26) 

1 = 1 

where Nj is an integer. The weights and scalar function are defined as: 

+1 = > (0 Nx +2 = II > +3 = SU P ||f “ RoEfRo I 

seR'CJT 

and V, = 1 ,V 2 =^ 1 ,^. 3 -|^/| 

y a is estimated as 

)-»>■(!") (27) 

i=l 

where N a = N\ +3, denote the estimate weights and ^n,(-) are some basis functions 
satisfying 


The compensating signal can now be defined as: 

Zv'n.^VnO') 

Z^n, Nl + V 


(28) 


(29) 


C.3.1. Stability Analysis 

For all control systems, and for adaptive control systems, stability is the primary requirement. In 
this section we use the Lyapunov’s stability theory to prove that the proposed new adaptive 
scheme in the above section is able to achieve the control objective. We choose a Lyapunov 
candidate function 
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( 30 ) 


V = V 1 +V 2 +V 3 


with 


r,=Vs, v 2 =-Ltr(fr, r fr,X 

2 2gj 


(31) 


^ 3 = 


2 Sn c „ 


< a i-c^&f{a>-c w J>) 


where gi > 0 and gn > 0 are free design parameters affecting weight learning rate and control 
performance, c m in > 0 is the minimum eigen value of REjR t . Its derivative along the vector field 
is 

V a =-kJ{R,E f Rj)s + s T W l T v l + s T 7 1 -s T R,E J Rju ca + V 2 + V 2 




Q +tr(w i T (y, 1 s T -g?lt) J 


with 


Q=s T r]-s T R n E f R 0 T u a , — a> 
£n 


(32) 


Before presenting the weight tuning algorithms, let us focus on Q. With u ca as in Eq. (30) and 
the upper norm bound on ij as in Eq. (25), we have 

C -Z«^ , 

Q*W Y - -s - —°> 


Z^M 


+v 


Sn 


= ( fi> - c ™^) r 


IX HIV* 


- O) 

Sn 


(33) 


+ VGJ 


Z^M +v 

/-i 

where co m = max (a),). Note that in deriving the above equation, we have used the facts that 

„ Z>//JMI 


a) 


Y < <o T y/ n 'Sn, and 


<1, Vv > 0 


Z *>11*11- 

;=1 

Combining Eq. (33) and Eq. (34), we come up with the following weight tuning algorithms: 
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( 34 ) 


W, = -<J,W, + g,y/,s 

S Vn, INI' Vn 

a = -cr u a + g„ - 


Consequently, 


5 >//, IHI +V 


V*-c min k 0 ||5|| 2 + -L tr (W 1 'W l ) 

g i 

+ —(a> - c mtn 6) T a + vco m 

8 ii 


Completing the square we can further rewrite Eq. (36) as: 


VZ-c mm k 0 \\s\?-^tr(W?W{) 
2 g, 

A\T 


5 —(fl> - c mi » r {co - C m ,» + 8 


2 §n c n 

= -A,V, - X 2 V 2 - X 3 V 3 + 8 


where 


8=va) m +-* l -t*W?W l )+-2- a) T a> 

2& 2gn 

and \ — 2c min & 0 , A^—g^ = cr n 

are constants independent of s. Note that Eq. (37) implies that ||s||, tr, and© are bounded. 
Furthermore, from Eq. (37), we have 

VZ-k^^f + S 

It is seen that v is negative as long as 


s\\ > 


8 N 

^0 C min J 


(35) 


(36) 


(37) 

(38) 


(39) 


(40) 


18 












Therefore, s is confined in the region n = 


. From previous lemma, we establish 



that the control error is bounded. 

C.3.2. Control Strategy with Alternate Gain 

An alternate control gain can also be used to achieve the above performance. This alternate gain 
enhances the performance of the control but the inverse is tedious to compute. The neuro- 
adaptive control for the system described as in Eq. (19) is given as: 

u = ^o K s ~ fm/ ~ w co) (41) 

where all the terms are as defined in the previous section. 


Properties 


1) The matrix £/ is a diagonal matrix that is positive definite. This fact can be verified by 
observing the range of e fi . 

2) Such matrix is nonsingular; therefore all the eigenvalues are non-zero. 

3) It is real symmetric; therefore all the eigenvalues are real numbers. The eigenvalues of 
such matrix can be given as : 


E f =diag(e fx . 

=> s fv ..£ fn are eigen values of E f 


(42) 


4) Matrix E is known from the literature as being invertible. As for matrix B, the control 
surfaces of the reusable launch vehicle are designed to control angular rate along each 
axis of the aircraft independently; therefore, the matrix B is invertible for all cases. 
Numerical studies for the aerodynamic model also show that B is always invertible. 

5) f Ro is full row rank, then there exist constants c min and cw* such that 


IHf s V 


5 < c m „ m 

max | | 


(43) 


where c miX > c min > 0 are the eigenvalues of the matrix R^E f R ^'. 


=min(e;g(£ / )) = f /mm >0 and = max («*(£,)) = £ fma > 0 (44) 


With the control of the above scheme, the original system Eq. (19) becomes 
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s=-^E f R;'s+W ^, -R^E f R 0 ' ] u ca +tj 


where W, =W,-W, is weight estimation error and 


7 = 5 + (7-R 0 £ / V)(/o + W i V i ) (45) 

is the resultant uncertainty to be attenuated by the second NN unit. Expressing rj in a similar 
fashion as done earlier, we have, |^| < y a 


N ,+2 

Ya = Z W > Y i (46) 

;=1 

where Nj is an integer. The weights and scalar function are defined as: 

^,+1 = d 0 , Wv I+ 2 = ||^;|| and 1^+, = 1, Y n>+2 = || y/, || 

y a is estimated as 

r.-£*i¥'iOl)“*Vi( F ) (47) 

i=l 

where N a = N\ +2, a>, denote the estimate weights and v/ n ,(.) are some basis functions satisfying 

V* ,(X,)^Yi (48) 

The compensating signal can now be defined as: 

Z^n^Vn I 1 ") 

U ca = "H*- * (49) 

2>n, Hl + V 

i -1 

C.3.3. Stability Analysis of the Control Strategy with Alternate Gain 

We choose a Lyapunov candidate function same as described in the previous section. Its 
derivative along the vector field is 

V a =-k 0 s T R i) E,R;'s + s T W. ,V, +s T n-s T R i) E / R;'u ca + V 2 + V 3 

*-*»*-. Mf + G 

with 
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( 50 ) 


Q = s T 7i-s T R 0 E f R 0 - ] u ca -—(w-c mm w) T w 

8 n 

Before presenting the weight tuning algorithms, let us focus on Q. With u ca as in Eq. (49) 
and the upper norm bound on 77 as in Eq. (45), we have 


--(w-c^w) w 


Sb 


= (H'-C m M’) 


Z^,H +v 


2>*H 


+ v 


1 i 

- w 

8n 


+ vw m 


(51) 


Note the expression for the signal Q is the same as Eq. (34). Thus, the rest of the stability 
analysis follows the previous section and we can prove that the control error is bounded in the 
region described as in Eq. (40). 

The control strategy change does not affect the performance of the controller in any way 
although the gain described in section A is easier to compute than the gain of control strategy 
described in the above section, but the latter method gives a better tracking precision. 


C.4. Simulation Results 

In this section, we present a series of simulations to illustrate the effectiveness of the proposed 
neuro-adaptive fault-tolerant control in tracking a desired angle of attack a, sideslip J3 and body- 
axis roll 0 of the UAV. 


Table 2: Numerical Values Used for Simulation 


Parameter 

Symbol 

Units 

Value 

Wing Area 

S 

m 

391.2 

Thrust 

T 

Lb 

27,000 

Mean Aerodynamic Chord 

c 

Ft 

11.32 

Roll Inertia 

l , 

kg.m 2 

1,997,922 

Pitch Inertia 


kg.m 2 

276,629,966.8 

Yaw Inertia 

h 

kg.m 2 

28,383,800 

UAV Mass 

M 

Kg 

103,756 

Wing Span 

B 

M 

29.2 

UAV Velocity 

V 

Mach 

0-15 
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The control strategy explained in earlier is utilized. Matlab is used for the simulation. The 
coefficients similar to the model X-33 are used. The simulation models include all non-linearity 
in the equations of motion, and the values of the parameters and coefficients are given as in 
Table 2 15,16,2 °. The payload mass is considered to be zero in the reentry phase and the mass of 
fuel is considered to be exponentially decreasing non-zero function with respect to time as 
m f = m f(ml j * e~',t- Time in seconds . 

We chose the basis functions for the first NN unit as: 


Vi(c u ,X,y) 


1 + e ^<M-cT ’ 


c;, = 0.2, 0.4, .... \\ x \\ = 


X = 0.5. 

The basis functions for the second NN unit are selected as 


^(y) = y,y=[W \\x\\ u . . l *],|*|« 


Fifteen NN basis functions were used for each of the NN units: y/eR 1 , j = 15. All the weights 

are initialized at zero in the simulation. The command values of a, (3 and <P are applied to the 

UAV are a =0.8 + 0.1 sin-, p' = 0.1 + 0.2e“', <f = 0.8 + 0.01 cos t . 

6 

The system faults for the numerical simulation are chosen to occur after 15 seconds and 
continue to appear till the end of the simulations. The faults represented in Eq. (14) with a = 0.5 
were chosen. 

The actuator faults occur after 30 seconds and also continue to appear all along till the end of 
simulations. These faults are modeled based on Eq. (16) with £fi chosen as 

fl, '^30 

Efi jl-J/*|sin (/-/)[,/> 30’ 

Figure 12 - Figure 18 shows the system response when the actuator power and sensors fail. 
Each figure has the actual and desired paths, the tracking error, and the action of the control, 
profile of faults and estimated weights for the first and second NN. It can be seen that the 
tracking is achieved within a small error. Note that the first NN weights remain unaffected by 
the fault occurrence. 
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Figure 12 - Actual and desired 
trajectories 
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Figure 13 - Control inputs 




Figure 14 - Trajectory Errors Figure 15- System Faults Profile 
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Figure 16 - Actuator faults profile 



Figure 17- Weights profile of first 
NN 



Time (Sec) 

Figure 18 - Weights profile of funneling NN 


The effect of number of neurons used for the neural networks plays a vital role in the 
effectiveness of the controller. It can be observed from the results presented above that the 
system reaches stability when 15 neurons were used. But the tracking precision can be improved 
by increasing the number of neurons in the NN. Figure 19 - Figure 23 shows the system 
responses when 25 neurons were used to represent a layer under system and actuator failures. It 
can be observed that even though the tracking precision is improved by a small amount, the 
controller exhibits chattering effect, which is a negative performance criterion for any control 
strategy. This trade-off is to be decided by the designer based on the nature of the plant. Design 
parameters like gj and gu aid in reducing this chattering but do not completely avoid. 
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Figure 19 - Actual and desired 
trajectories 
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Figure 20 - Trajectory errors 
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Figure 21- Controls inputs 


Figure 22 - Weights profile of first NN 
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Figure 23 - Weights profile of the funneling NN 
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Appendix I -Nomenclature 


b 

c 

m 


<7 

S 

T 

V 

F 


L, M N 


S'faA 

P 

c 


= UAV wingspan 

= Wing Mean Reference Chord 

= UAV Variable mass 

= Roll, Pitch and Yaw rates/body fixed frame 

= Free-Stream dynamic pressure 

= Wing Reference Area 

= Thrust 

= UAV variable velocity 
= Aerodynamic Force about the body frame 
= Aerodynamic Rolling, Pitching and Yawing Moments 
= Angle of Attack and Side-slip 
= Roll, Pitch and Yaw Euler angles 
= Elevator, Aileron and Rudder deflections 
= Air Density 

= Aerodynamic Force Coefficient 
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